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The hair cells of the vertebrate inner ear convert mechanical stimuli to electrical signals. Two 
adaptation mechanisms are known to modify the ionic current flowing through the transduction 
channels of the hair bundles: a rapid process involves Ca 2+ ions binding to the channels; and a 
slower adaptation is associated with the movement of myosin motors. We present a mathematical 
model of the hair cell which demonstrates that the combination of these two mechanisms can 
produce 'self-tuned critical oscillations', i.e. maintain the hair bundle at the threshold of an 
oscillatory instability. The characteristic frequency depends on the geometry of the bundle and 
on the Ca 2 + dynamics, but is independent of channel kinetics. Poised on the verge of vibrating, 
the hair bundle acts as an active amplifier. However, if the hair cell is sufficiently perturbed, other 
dynamical regimes can occur. These include slow relaxation oscillations which resemble the hair 
bundle motion observed in some experimental preparations. 
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Introduction 

Hair cells of the inner ear detect mechanical stimuli 
by deflections of the hair bundle, which open tension- 
gated transduction channels in the cell membrane to ad- 
mit cations from the endolymph. The current is mostly 
composed of potassium io ns, but also inc l udes a small 
quantity of calcium ions l)Lumpkin et all Il997|) . Ex- 
periments in which hair bundles are held at a fixed de- 
flection, using a microneedle, have revealed two adap- 
tation processes that modify th e transduction current , 
causing it t o dec li ne with time pfolt and Coreyl 
[ Holt et all 120021 iHoward and Hudsoethl Il987t 
IWu et all 1 1 999(l . The first, which can occur on a sub- 
millisecond time scale, is dependent on the concentration 
of Ca 2+ in the endolymph. It is believed to be a conse- 
quence of Ca 2+ entering the hair bundle and binding to 
the transduction channels, making them more likely to 
close. The second, which takes up to hundreds of mil- 
liseconds, is also dependent on Ca 2+ concentration, and 
is thought to be caused by the movement of adaptation 
motors attached to the transduction channels, which ad- 
justs the tension that gates them. 

A number of experiments have demonstrated that 
hair bundles can generate active oscillations, both in 
preparations of dissected tissue llBenser et all 19931 
ICrawford and Fettinlace]. Il985t iMartin et all l2000|) and 
in vivo llManTeT^talT " 2001^1 . It has been shown on 
theoretical grounds that the active oscillatory mecha- 
nism can provide an explanation for the high sensitiv- 
ity and sharp frequen c y selectivity of auditory receptors 
iCamalet et all l2000t IChoe et all Il998t lEguiluz et all 
120001 lOsneck et alll20oiF if the system is poised on the 
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verge of spontaneous oscillation. This critical point cor- 
responds to a Hopf bifurcation of the active dynamical 
system. It has been suggested that a feedback mech- 
anism operates to maintain the system at the critical 
point where it is most sensitive, leading to the concep t 
of self-tuned critical oscillations. l|Camalet et all [2000). 

The purpose of this article is to provide a coherent 
model for the active amplifier in hair cells, based on the 
combined action of the two known adaptation mecha- 
nisms. In our model the interaction of Ca 2+ ions with 
the transduction channels generates an oscillatory insta- 
bility of the hair bundle. Working on a slower time scale, 
the adaptation motors act to tune this dynamical system 
to the critical point, at which it just begins to oscillate. 
The two adaptation mechanisms together constitute a 
dyna mical system display ing a self-tuned Hopf bifurca- 
tion l)Camalet et allfeoOOj) . Poised at the oscillatory in- 
stability, the hair bundle is especially responsive to faint 
sounds. 

Our model differs from two published models of hair- 
bundl e oscillations. Choe, Magnasco and Hudspeth 
(1998) have previously proposed that channel reclosure 
induced by Ca 2+ ions can generate a Hopf instability. 
In their scheme the oscillation frequency is controlled by 
the kinetics of the transduction channel. In our model 
the frequency of spontaneous oscillations is governed, in- 
stead, by the architecture of the bundle and by the rate 
at which the calcium concentration equilibrates. This is 
more consistent with the finding that the rate at which 
the transduction channel switches is faster than the char- 
acteristic frequency, especial l y in low-frequency hai r cells 
l|Corev and Hudspethl . ll983]) . ICamalet et al 1 il2000l) have 
suggested that the Hopf instability is generated by mo- 
tor proteins in the kinocilium, and that regulation to 
the critical point is accomplished by a feedback mech- 
anism involving calcium. Their model is thus the con- 
verse of that presented here, but shares the underlying 
feature of self-tuned criticality. The current model ac- 
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cords better with the observation that some hair bundles 
lack a kinocilium and others do not change their behav- 
ior significant l y whe n the kinocilium is interfered with 
l|Martin et alJ.l2003|) . 

In addition to demonstrating that the two adaptation 
processes can combine to provide an active amplifier, we 
derive expressions that indicate how properties like the 
characteristic frequency of critical oscillations, and the 
sensitivity of response to sound waves, depend on param- 
eters like the bundle height and the number of stereocilia. 
We also investigate how the dynamics can be modified if 
conditions are perturbed, and show that the bundle can 
exhibit slow relaxation oscillations which differ in many 
aspects from critical oscillations, but resemble the active 
movements recently obs erved in experiments on hair cells 
from the frog sacculus ijMartin and Hudspethl . 11999). 

Model 
Hair bundle 




FIG. 1 (a) Schematic representation of a hair bundle, (b) 
Three-state model of the transduction channel. Note that 
Ca 2+ binds to the channel when it is already in a closed state. 
The presence of Ca 2+ therefore does not directly cause chan- 
nel re-closure, but it does stabilize the closed state and delay 
re-opening. 

Hair bundles in different species, or from hair cells with 
different characteristic frequencies within the same organ 
of thte internal ear of a particular animal, can vary in 
shape and size, and also in the way t hat they are cou - 
pled to external mechanical vibrations l|Hudspetbl .ll997'). 
Generally, one can distinguish between free-standing hair 



bundles (e.g. some reptiles, mammalian inner hair cells) 
and those who tips are embedded in an accessory struc- 
ture (e.g. other reptiles, birds, amphibial saccular hair 
cells). In experiments on dissected tissue, the accessory 
structure is usually removed. In this paper we will re- 
strict our discussion to the simplest situation of a free- 
standing hair bundle. A bundle that is directly connected 
to the tectorial membrane could be described in a similar 
fashion, but using modified expressions for the stiffness, 
the friction and the coupling to the external signal. 

The hair bundle (see Fig. ^,) is composed of N stere- 
ocilia which slope up against each other, the longest of 
which has length L. Typically, N an d L are inversely 
relate d, with the range N = 20 — 300 l|Holt and Corevl 
2000) corresponding to L = 30—2 /im. Each stereocilium 
is joined to the next by an elastic filament - the tip link 
which connects to a t ension-gated transduction ch annel 
in the cell membrane ( Corcv and Hudspethl . Il983t for a 
review see iMarkin and Hudspethl II995I) . The tension in 
the tip link can be adjusted by a group of myosin mo- 
tors, which are connected to the transduction channel 
and mo ve it along the actin filaments inside the stere- 
ocilium l)Hudspeth and Gillespiei Il994jh 

To set up a coordinate system, we consider first a free- 
standing bundle. Within a stereocilium, we label the lo- 
cation at which the stereocilium touches its shorter neigh- 
bor O. The position of the channel (and the attached 
motors) relative to O is denoted by Iq + xm, where Iq 
is the length of a tip link at rest and % its extension 
under tension. When the tip of the bundle is deflected 
through a lateral displacement X, each stereocilium piv- 
ots at its base and adjacent stereocilia are sheared by a 
displacement 'yX. The geometric factor 7 may be written 
as 7 w b/L, where b ss I fim is t he distance between the 
roots of neighboring stere ocilia ijHoward and Hudspethl . 
lf988tlHoward et allll988|) . We label the new location at 
which the stereocilium touches its shorter neighbor O'. 
The position of the channel relative to O' is denoted by 
lo + x. It follows that 

jX = x - xu ■ (1) 

The movement of the bundle is countered by an elastic 
restoring force, which has two contributions. First, the 
deformation of the stereocilia at their base provides an 
effective bundle stiffness Kgp = kN/L 2 , where k is the 
pivotal stiffness of a ste reocilium. The measured value 
k ft 1.5 x I0~ 16 Nm/ra d l|Crawford and Fettiplacelll985l 
iHudspeth et all |2£)00) indicates that Kgp varies from 
O.OIpN/nm for tall bundles to fOpN/nm for short ones. 
The second contribution, which scales in the same way 
with bundle size and is of similar magnitude, comes from 
the tension Ttl in the tip links. Movement of the bundle 
is also resisted by the viscous drag of the surrounding 
fluid which leads to a friction coefficient £ ~ rjL for tall 
bundles, where i] is the viscosity of the endolymph. Thus, 
when an external force F is applied to the tip of the bun- 
dle, its equation of motion is 

(X = -N-yT TL (x, C) - K SP X + F , (2) 
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(throughout the paper we use the convention X(t) = 
^-). Here, we have written T TL = T TL (x, C), antici- 
pating that the tension in the tip links depends on the 
local intracellular Ca 2+ concentration C, as well as vary- 
ing with the position x. For future reference, we note 
that the mechanical relaxation time of the hair bundle, 
Tmcch = (/Ksp, varies from a few microseconds to tens 
of milliseconds. 



Transduction channel 

The mechanically-gated transduction channel is the 
central component of the model, since we postulate that 
hair bundle oscillations are driven by the force that is 
generated as the channels change conformation stochasti - 
cally, as first suggested bv lHoward and Hudspethl fl988). 
It is generally assumed that transmembrane channel pro- 
teins can exist in a number of discrete states, some of 
which are closed and some open. As indicated in Fig.^D, 
we additionally suppose that the channel protein incor- 
porates a lever arm which amplifies the small structural 
changes that occur when the channel switches state, and 
thereby appr eciably modifi e s the tension in the adjoin- 
ing tip link ijMartin et all l2000|) . If pi is the proba- 
bility that the channel is in state i, and di is the po- 
sition of the lever arm in that state, then the aver- 
age tension may be written Ttl = Xh-^tlOe — di)Vi- 
Here, -RTtl is the elastic constant of the tip link (or 
possibl y that of a more compliant element in series 
with it ijKachar et alll2000^ ). whose value obtained from 
micromechani c al me asurements is Ktl ~ 0.5pN/nm 
((Martin et all l2000(l . In the steady state the proba- 
bilities pi = Wi/ (y2j Wj J are given by the Boltzmann 

factors Wi — exp [— (iif TL (j; — di) 2 — G®)/k B T], where 
G° is the free energy of state i. Measurements of the 
transduction current as a function of bundle displace- 
ments in the frog's sacculus are well fitted with a three - 
state model of the channel IjCorev and Hudspethl . Il983). 
with two closed states (i = C1,C2) and one open state 
(i = O). The fit indicates that the lever-arm move- 
ment is about 5 times greater for the transition CI — > C2 
than for the transiti on C2 — * O, while more recent data 
IIMartin et aH|2000) suggests that the total movement is 
about 8 nm. We therefore assign the lever arm displace- 
ments dci — Onm, dc2 = 7nm and do = 8.5 nm. The 
free energies G° depend o n the local Ca 2+ concentratio n 
C within the stereocilia ijHoward and Hudspethl . Fl988'l. 
which is strongly infl uenced by the entry of Ca 2+ ions 
through the channel ijLumnkin and Hudspethl . Il995ll . A 
higher Ca 2+ concentration favors closure of the channel, 
thereby providing a negative feedback. We model this 
depende nce empirically, according to experimental obser- 
vations ijCorev and Hudspethl . 11983^1 . as an effective shift 
in the channel gating point proportional to the logarithm 



of the Ca concentration 

Pi(x,C) = Pi (x-x s (C)) , x s (C)=Dln — , (3) 

where the constant D has an em pirical value of about 
4 nm in the bullfrog's sacculus l(Corev and Hudspethl . 
1983) and Co is the Ca 2+ concentration in the free- 
standing bundle. One possible interpretation of this func- 
tional form, indicated in Fig. is that Ca 2+ binds to 
state C2, causing it to switch to state CI and thereby sta- 
bilizing the closed state. We emphasize that our model 
does not rely on the channel having three states. A two- 
state channel (eliminating state C2) would give substan- 
tially similar results. 

An important aspect of our model is that we as- 
sume that the channel kinetics are rapid compared 
to the bundle oscillations and the calcium dynamics. 
The evidence for this assumption comes, for instance, 
from saccular hair cells that detect frequencies up to 
150 H z and have channel op ening times around 100 fis 
l|Corev and Hudspethl . Il983|i . We can therefore consider 
the different channel states to be thermally equilibrated 
with one another, and chemically equilibrated with the 
Ca 2+ ions. Then we can split the tip link tension Ttl into 
two components: K^x, the passive elastic force that the 
tip link would contribute in the absence of the transduc- 
tion channel; and /, the active contribution arising from 
the channel switching states, also called gating force: 

T Th (x,C)=K Th x + f(x-x s {C)), (4) 

where f(x,C) = —J2i^ThdiPi(x,C). A relation- 
ship between the gating force and the channel open- 
ing probability, which holds in any model with a single 
open state, reads fix C) = fe B r gtofeg]3 _ K TL d 
llvan Netten and KrosL 12000ft . 



Calcium dynamics 

The dynamics of the calcium concentration in a stere- 
ocilium is influenced by a number of processes. As 
well as the influx of Ca 2+ through the transduction 
channel, there is an outflux due to calcium pumps. 
The ions also diffuse freely, or in association with mo- 
bile buffers, and can additionally bind to fixed buffers. 
A detailed model of th i s dyn amics was developed by 
iLumpkin and Hudspethl l|l998f) . but for our purposes it 
suffices to use a simplified linearized model in which the 
Ca 2+ concentration on the inner side of the channel obeys 
the equation 

C = -X[C-C B -C M p {x,C)] . (5) 

Here Cb is the steady-state concentration which would 
occur if the channel were blocked, originating from dif- 
fusion from other parts of the cell; Cm is the maxi- 
mal additional contribution from influx through an open 
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transduction channel, and is of order 10 fj,M. The re- 
laxation rate of calcium concentration fluctuations, A, is 
an important parameter in our model. We expect the 
relaxation time r c h e m = 1/ A to be of the order of mil- 
liseco nds, corresponding to the m easured fast adaptation 
time (fRicci and Fettinlacei 119971) . Note that the slow 
Ca 2+ kinetics (along with fast channel swi tching) is a 
key di fference between our model and that of lChoe et al.l 
(1998), where the converse situation is assumed. The 
slow Ca 2+ dynamics might seem at odds with the fact 
that the binding sites on the channel are reached by the 
incoming ions very quickly. But in our model, the chan- 
nel is in state C2 when Ca 2+ binds, i.e. it is closed 
(Fig. It therefore typically senses a concentration 

which is determined by diffusion within a larger volume, 
and which also depends on the interaction of Ca 2+ with 
fixed and mobile buffers. Becau se a Ca 2+ ion binds to one 
of th e buffers in less than 1 fis l|Lumpkin and Hu dsoc th, 
1998), the concentration of free Ca 2+ equilibrates with 
that bound to the buffers within a few microseconds 
of channel closure O^C2. The time-scale of Ca 2+ re- 
laxation is therefore determined by buffering, diffusion 
to distant parts of the stereocilium and the activity of 
Ca 2+ pumps. Indeed, the influence of buffers on the 
fast adaptation process has been observed experimentally 
l)Ricci and Fettiplacelll997h . 



Adaptation motors 

Experiments i n which a vertebrate hair bundle is sud- 
denly displaced ((Howard and Hudspethl . lT987llWu et, all 
1999) show a slow adaptation of the channel current, 
on the scale of tens of milliseconds. The transduc- 
tion channel is attached to the actin cortex within a 
stereocilium by a number of myosin 1C (formerly 1/3) 
motors. There is strong evidence t hat these molecu- 
lar motors mediate slow adapatation IIAssad and Corevl 
119921 iGillespie and Corevl Il997t llloit et all l2002t) . pre 
sum ably by maintaining the prope r tenstion in the tip 
link ([Hudspeth and GillespieL 119941) . For simplicity, we 
assume that the motors have a linear force-velocity rela- 
tion, so that 



xm = v M 



1 



M 



(6) 



where tim is the number of motors bound to the actin, 
Fm is the stall force of an individual motor and vm is 
its zero-load velocity. We take values typical for a motor 
protein, Fm ~ 1 pN and % 0.3 /an/s and suppose that 
there are Nm ~ 20 m otors per channel (experiment al 
data suggest Nm < 130 l|Gillespie and Hudspethl .LI993)). 

Calcium has been shown to influence the action of 
the a daptation motor ijHacohen et all Il989l iRicci et ail 
1998), possibly by increasing its rate of detachment from 
actin. The number of bound motors is therefore not con- 
stant, but obeys the dynamical equation 



where u>m ~ 200 s _1 is the myosin binding rate and g(C) 
is a function describing the dependence of the detachment 
rate on the Ca 2+ concentration. 



Results and Discussion 

Linear stability analysis 

In order to examine the dynamical stability of the hair 
bundle in the absence of a driving force, F = 0, we make 
use of the fact that the adaptation motors operate on 
a slower time scale than both the mechanical relaxation 
of the bundle and the relaxation of the calcium concen- 
tration. We can therefore neglect their motion for now, 
and suppose that Xm is fixed. This leaves two coupled 
first-order differential equations of motion: Eq. [S] (sup- 
plemented with Eq.[3J) for the rate of change of Ca 2+ con- 
centration; and by combining Eqs.^|2]&0J the following 
equation for the motion of the transduction channels: 

x = - J{iV 7 2 [Ktlx + f(x- x s {C))\ + K SP (x - x M )} ■ 

The system has one or more fixed points, given by 
x = and (7=0, whose stability can be determined 
from the Jacobian matrix 




n M = lo m \-n M g{C) + (N M - n M )} 



(7) 



(9) 



which has trace and determinant 

Tr J = -i [K B + N 7 2 f] — A (1 + /z) , (10) 

det J - ~ [X B (l + A*) + ^7 2 /'] , (11) 

(throughout the paper we use the convention f'(x) = j^, 

p' (x) = ^ andx' s = Here, K B = K sp +N 7 2 Ktl 

is the combined passive elasticity of the bundle due to 
the stereocilia and the tip links; and \x — Cmx' s p'q is 
a dimensionless combination of variables which has a 
value close to unity for the choice of channel parame- 
ters specified above (for Cb — it can also be expressed 
as (i = D(f + KThdo)/kBT). A fixed point is stable if 
both eigenvalues of J, evaluated at the fixed point, hav e a 
negative real part. This is equivalent l)Strogatzlll994^ to 
the two conditions: det J > 0, and Tr J < 0. The number 
of fixed points is most readily appreciated by consider- 
ing the form of Eq. [H] for the velocity function of 
displacement x (see Fig. Two distinct situations can 
arise depending on the strength of the calcium feedback: 
(i) Consider first the case of weak feedback, so that 
the Ca 2+ concentration C in the stereocilia is approxi- 
mately constant. Typically x(x,C = const) varies non- 
monotonically with x. This is a consequence of a region 
of effective negative elasticity, caused by the redistribu- 
tion of the t ransduction channel states when the bun- 
dle is moved ijHoward and Hudspethlll98^:lMartin et all 
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FIG. 2 Velocity x as a function of displacement x for a fixed 
location of the motors xm- If the Ca 2+ concentration is held 
constant, x(x,C = const) typically has a region of negative 
slope (solid lines, for three different hypothetical values of 
C). For a certain range of values of C, the system can then 
have three fixed points. But if the Ca 2+ concentration is 
allowed to adjust, the curve x(x,C = C(x)), where C(x) is 
the steady-state solution of Eq. [H] rises monotonically with x 
(dotted line). There is then a single fixed point. The channel 
parameters are G C1 = 0, Gc 2 = 60 zj, G% = 70 zj, dci = 
Onm, dc2 = 7nm, do = 8.5 nm and D — 4nm. With these 
channel parameters KtlcIq w 9fcsT (in a simplified two-state 
channel model neglecting the stereociliary stiffness Ksp the 
condition for the occurrence of a negative slope is KtlcIq > 
ik B T). 



2000). At large displacements - both positive and neg- 
ative - the tension Ttl in the tip links is determined 
by their passive elasticity. But at intermediate displace- 
ments, the occupancy of the closed state CI declines 
rapidly with increasing x; the resulting force /, caused 
by the swing of the lever arm, is negative and reduces the 
tension in the tip links so that Ttl decreases with increas- 
in g x. This effect has been term ed 'gating compliance' 
bv lHoward and Hudsnethl ()l987j) . As indicated in Fig.^J 
there will then be either one or three fixed p oints. Ac- 
cording to the index theorem l)Strogatzt Il99l Sect. 6.8) 
out of 2n+ 1 fixed points, n will be saddle-points. There- 
fore, if the system has three fixed points, the central fixed 
point is a saddle point and the other two may be either 
stable or unstable. The number of fixed points depends 
on the values of the channel parameters, but also on the 
motor position xm- Their stability depends additionally 
on the values of the dynamical model parameters (£ and 
A). A typical example of the location and nature of the 
fixed points, as a function of the motor position xm, is 
shown in Fig.UJi; in this particular case, the bundle dis- 
plays a region of bistability. We shall return to this situ- 
ation later, to discuss the dynamical behavior of the hair 
bundle in recent in vitro experiments. 

(ii) When the calcium feedback is operational, to de- 
termine the fixed points it is appropriate to consider the 
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FIG. 3 Fixed points of the displacement 1 as a function of 
the location of the adaptation motors xm. Trajectories in 
the plane (a;, C) are indicated schematically below, (a) The 
situation with a weak Ca 2+ feedback. The system displays a 
variety of dynamical regimes, including a region of bistability. 
(b) With strong Ca 2+ feedback there is a single fixed point 
which becomes unstable (and encircled by a stable limit cycle) 
for an intermediate range of values of im. We expect the hair 
cells to use this regime in vivo. 



curve x(x,C = C(x)), where C(x) is the steady-state 
Ca 2+ concentration at fixed x, given by C = in Eq. [5] 
If the feedback is sufficiently strong, x varies monoton- 
ically with x (Fig. |2J ; the readjustment of the channel 
states caused by the binding of Ca 2+ ions provides a pos- 
itive contribution to the active force /, eliminating the 
region of negative elasticity. The condition for this to 
occur is that det J > along the curve C{x), or 



ft > 



NY 



/'-I 



(12) 



In this case, there is a single fixed point (Fig. [3b), 
which may be either stable (if Tr J < 0) or unstable (if 
TrJ > 0). In the latter situation, the bundle under- 
goes limit-cycle oscillations. Significantly, the stability 
depends on the location of the motors xm, in addition to 
the values of the other model parameters. Typically the 
critical point at which Tr J = is a supercritical Hopf 
bifurcation, which means that the system is continuously 
controllable, passing smoothly from the quiescent to the 
vibrating state as indicated in Fig.|Ij Thus by moving 
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FIG. 4 Displacement i as a function of the location of the 
adaptation motors xm, for the proposed in vivo situation cor- 
responding to Fig. Eh. Between the two bifurcation points, x 
oscillates between a minimum and a maximum value shown 
by the two curves. The motor-mediated self-tuning mecha- 
nism (see text) maintains the system at a working point close 
to the first Hopf bifurcation. This result was obtained from 
a stochastic simulation with full channel kinetics, but in the 
absence of external (thermal) noise. The upper and lower line 
show the average maximum and minimum values of the dis- 
placement during one period of oscillation. Note that even in 
the quiescent phase the amplitude is non-zero, due to channel 
noise. The noisy motion in the quiescent phase has a very 
different spectrum from the oscillating phase, however. 

along the stereocilia, the motors can determine whether 
the bundle stays still, or with what amplitude it oscillates 
spontaneously. 

Frequency of limit-cycle oscillations 

We propose that the regime of strong calcium feed- 
back, in which the bundle can oscillate spontaneously, is 
the regime that is relevant for auditory hair cells in nor- 
mal physiological condititions. The control parameter 
that determines the stability of the bundle is e — | Tr J 
and the Hopf bifurcation occurs at e = 0. It is worth em- 
phasizing that within the scope of a non-inertial model 
with two variables (x and C), in order for the bundle 
to oscillate spontaneously, it is essential that either the 
bundle displays a region of negative elasticity at constant 
Ca 2+ concentration, or that there is a positive feedback 
in the Ca 2+ dynamics. The second option would imply 
that an increase in the intracellular Ca 2+ concentration 
would enhance the Ca 2+ inflow, which contradicts the 
present experimental evidence. Our model only allows 
for the first option and from Eq. it can be seen that 
the critical condition e = requires that Ks+N^f 2 /' < 0, 
i.e. the passive bundle stiffness must be more than com- 
pensated by the gating compliance. Negative elasticity 
has been mea s ured in hair cells of the bullfrog sacculus 
ijMartin et alll2000(l . and has been interpreted as being 



due to gating compliance ijHoward and Hudspethl . Il987j) 
as discussed above. 

For small, positive values of e, the bundle ex- 
ecutes small-amplitude limit-cycle oscillations of fre- 
quency log = vdetJ. From Eq. ^2 it follows that the 
characteristic frequency of the hair bundle scales as 



Uc ^ y^A^p. oc y'A^ (13) 

Note that the period of oscillation Tc is approximately 
given by the geometric mean of the two relaxation times, 
Tc ~ ^/r mcc hT c hem- Little is known about the variability 
of the calcium relaxation rate A. If it is constant from cell 
to cell, the characteristic frequency is determined only by 
the bundle geometry. For the range of bundle sizes found 
e.g. in the chick's cochlea (L = 1.5 — 5.5/im, N = 300—50 
l|Tilnev and SaundersL ll983)). the frequency of the short- 
est bundles would then be approximately 17 times greater 
than the frequency of the tallest bundles. This is some- 
what lower th an the measured fre quency range (about 
50 Hz to 5 kHz l)Manlev et allll99irj ). We note, however, 
that to ensure the existence of an oscillatory instability 
for all hair cells, e = from Eq. 1101 must be fulfilled for 
each cell. In addition, a strong mismatch between r moc h 
and Tchom would seriously limit the range of amplitudes 
for which active amplification can be maintained. This 
suggests that A may vary with the geometry of the cell; 
ideally it would scale as A ~ Kb/(, i-e. T c h cm ~ Tmcch- 
Such variation could be achieved by a differential level of 
expression of calcium pumps or different buffer concen- 
trations. In this case, the range would be enhanced to 
cover frequencies differing by a factor of 200. In what fol- 
lows, we shall assume that this scaling (which is equiva- 
lent to A ~ ll>c) holds. Note that in contrast to the m odel 
proposed by Choe, Magnasco and Hudspeth l)l998j) no 
variation of transduction channel kinetics is necessary in 
order to achieve a broad range of characteristic frequen- 
cies. Instead a simultaneous variation of the hair bundle 
geometry and the Ca 2+ relaxation rate is required. 



Critical oscillations 

It has been proposed on general theoretical grounds 
that the frequency selectivity of hair bundles, as well as 
their sensitivity to weak signals, is con ferred by their 
proximity to a dynamica l instability l|Camalet et all 
l2000t lEguiluz et all l2000|) . Poised on the verge of os- 
cillating at a characteristic frequency ljc, a hair bundle 
is especially responsive to periodic disturbances at that 
frequency. To ensure robust operation, a control sys- 
tem is required which maintains the hair bundle close 
to the critical point at which it just begins to vibrate. 
The general concept of a self-tuned Hopf b ifurcation was 
therefore introduced ijCamalet et all l2000() , whereby an 
appropriate feedback mechanism automatically tunes the 
system to the vicinity of the critical point. For the mo- 
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ment, we will assume that this mechanism works effec- 
tively, so that each hair cell is poised at a working point 
very close to its critical point, as indicated in Fig.0] 

In the vicinity of the bifurcation the limit cycle can be 
described by expanding the equations of motion about 
the fixed point as a power series, keeping terms up to 
third order. It is convenient to introduce a complex vari- 
able 



Z = Ax + a AC, 



(14) 



where Ax and AC measure the difference of x and C 
from their respective values at the critical point. With 
the choice 



J. 



22 - 6 + HOC 
Jl2 



(15) 



For sound frequencies close to the characteristic fre- 
quency of the bundle, u « loq, and amplitudes larger 
than the thermal noise, the response of the system is 



ReB 



(19) 



where we have assumed that the bundle is poised pre- 
cisely at the Hopf bifurcation, e = 0. 

Using the estimate Re B ~ -wc / d 2 ^ (see Appendix) we 
obtain 



Ax . 



^vAd 2 



vAd 2 Q b/L 



(20) 



where J i2 = XCmv'o an< ^ ^22 = — A(l + fjt) (Eq. 0, the 
equations of motion can be written as 

Z = {e + iuj G )Z + Q + B\Z\ 2 Z + C , (16) 

where Q represents all quadratic terms and C all remain- 
ing cubic terms. The expression for the coefficient B is 
given in the Appendix; for typical parameters, ReB is 
negative, in which case the bifurcation is supercritical 
and therefore continuously controllable. To leading or- 
der, the amplitude of the variable Z is then given by 
A = yf—e/ReB, while all other terms cau se only pertur- 
bations to the trajectory l)Wigginsl 119901. The solution 
for spontaneous critical oscillations therefore reads 

Ax( f ) = ^i, AC(t) = l -^l, Z(t)=Ae^ 



Im a 



Im a 



(17) 



and the oscillation amplitude of Ax is f^-A. 



This is the strongly compressive response characteristic 
of a H opf resonance l)Camalet et all l2000t lEguiluz et all 
2000). Displacements of the transduction channel which 
produce a resolvable response, Ax = 1 — 10 nm, typically 
correspond to wave amplitudes A = 1 — 1000 nm in the 
fluid. 

For weak stimuli, it is important to take the effects 
of noise into consideration. The buffeting of the hair- 
bundle by the Brownian motion in the surrounding fluid 
will render the spontaneous critical oscillations noisy. In 
this regime, signals that yield a response that is smaller 
than the spontaneous motion can nevertheless be de- 
tected, owing to pha se-locking of the noisy oscillations 
ijCamalet et allfeOQOj) . The ability to detect phase lock- 
ing in neural signals despite no increase in the firing rate 
has been experimentally co nfirmed in a number of species 
l)Narins and Wagnerlll989]) . In the Appendix we provide 
an argument for the threshold of detection which suggests 
that for high frequency bundles containing many stere- 
ocilia, the smallest wave amplitude that can be detected 
may be as small as A sa 0.1 nm. 



Critical response to sound waves 

Sound waves or other vibrations entering the inner ear 
cause flow in the fluid surrounding the hair bundles. If 
the (unperturbed) fluid motion at the height of the bun- 
dle tip has amplitude A and frequency d, the force ex- 
erted on the bundle tip in Eq. [21 is F = Auj(tiow cos(wt). 
Here, Cflow ~ i]L is the friction coefficient of a stationary 
bundle in a flowing fluid. Note that £flow is not pre- 
cisely equal to £, which is the friction coefficient of the 
bundle moving in a stationary flow, because the latter 
also includes a contribution from shear flow between the 
stereocilia. This contribution scales as ~ r\L x N^ 2 and 
is negligible for tall bundles but might become significant 
for short thick bundles. To preserve generality we write 
Cflow/C = v i where v is a constant of order unity which 
might vary with bundle geometry. 

The external force leads to an additional driving term 
on the right hand side of Eq. 1161 given by 

T = ^A^ue iut . (18) 



Self-tuning mechanism 

In order to achieve the combination of frequency se- 
lectivity, amplification of weak signals and wide dynamic 
range, the system has to be tuned to the close proxim- 
ity of the Hopf bifurcation (Tr J = in Eq. ITDf> . We 
propose that this is accomplished by a feedback control 
mechanism which is mediated by the adaptation motors. 
The basic requirement is that when the hair bundle is 
still, the force exerted by the motors gradually increases 
so that they advance up the stereocilia (increasing %), 
pushing the system into the unstable regime (see Fig.0J; 
but when the bundle is oscillating, the motor force grad- 
ually declines, so that the motors slip back down the 
stereocilia, returning the system to the quiescent state. 
It would appear that some type of high-pass or band- 
pass filter, which can detect whether or not the bundle 
is oscillating, must form an integral part of the feedback 
mechanism that controls the motors. 

It has long been known that the hair cells of some ver- 
tebrates have, in their cell membranes, a system of potas- 
sium channels (BK) and voltage-gated calcium channels 
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that display a resonant response t o oscillating input cur- 
rents ijHudspeth and LewisUl988|) . It has been suggested 
that their purpose is to sharpe n the tuning of the signal 
passe d to the auditory nerve ijCrawford and Fettiplacet 
1980). It has alternatively been postulated that this sys- 
tem may itsel f gener ate self-sustained Hopf oscillations 
i|Ospeck et all l200l|) and could be the basis of active 
bundle movement. We put forward a further possibil- 
ity: One role of this system of channels is to detect bun- 
dle oscillations and to control the adaptation motors via 
a flux of calcium ions into the hair bundle (in fact, it 
is more likely that the motors are regulated by Ca 2+ - 
controlled secondary messengers diffusing into the hair 
bundle, rather than by Ca 2+ directly, but this would not 
qualitatively alter the following argument). 

We model the filtering characteristic of the electrical 
oscillator by a Lorentzian transmission function 



V((jj) = Vo + /input M-Rf- 



(21) 



where V is the receptor potential, Rp is a channel resis- 
tance and /input oc po is the transduction current. The 
mechanism operates best if the filter frequency wf equals 
the characteristic frequency loq of hair-bundle oscilla- 
tions, but this matching has to be encoded genetically 
or achieved by some mechanism beyond the scope of our 
model. With the value Q s» 10 measured experimentall y 
in the bullfrog's sacculus IjHudspeth and Lewisl [1988J), 
the filter distinguishes sufficiently well between the os- 
cillating component of the transduction current and the 
underlying background level. The voltage-gated calcium 
channels then provide a Ca 2+ flux varying (for exam- 
ple) exponentially with the membrane potential, which 
in turn determines the offset concentration in the stere- 
ocilia Cb , which appears in Eq. [5j We model it with the 
following differential equation 



Cb = lea 



1 



l + exp(-g(F ~V Ca )/k B T) 



(22) 



where the parameters q and Vc a describe the voltage sen- 
sitivity of the Kca channels, /c a is a measure of these 
channels' contribution to the Ca 2+ flux in the region close 
to the transduction channel, and Ab is the relaxation rate 
of the Ca 2+ current. For the operation of tuning mech- 
anism it is essential that Ab *C loc- The contribution 
of Cb is insignificant when the bundle is quiescent, but 
grows rapidly if the bundle starts to oscillate. 

The Ca 2+ flowing into the bundle from the cell body 
regulates the oscillations in two ways. First, it reduces 
the control parameter of the bifurcation directly, Eq. 1101 
Second, Ca 2+ ions binding to the myosin motors enhance 
their rate of detachment, via the function g(C) in Eq. 
thereby regulating the total force that the team of motors 
produces. Here we take g(C) oc C 3 ; the precise functional 
form is not very important. 

As shown in Fig. this feedback control mechanism 
can successfully tune the system to the immediate vicin- 
ity of the Hopf bifurcation. In the absence of a sound 




FIG. 5 Self-tuned critical oscillations. In the absence of a 
stimulus, the bundle performs noisy oscillations with a small 
amplitude. A strong external stimulus elicits a response that 
can detune the system for some time. Parameters values: 
N = 50, A" TL = 0.5pN/nm, K S p = 0.15pN/nm, £ = Cflow = 
0.65 pN ms/nm, 7 = 0.14 and A = 40 nm ( the values of other 
parameters are listed in the last section of the Appendix ). 



stimulus, a bundle executes noisy, self-tuned critical os- 
cillations of low amplitude. If a strong stimulus is then 
applied, the bundle displays the characteristic response 
given by Eq.[2U] The significant flow of Ca 2+ into the cell 
elicited by this response detunes the system, so that when 
the stimulus is cut, the bundle falls still. Gradually, how- 
ever, the motor-mediated self-tuning mechanism returns 
the system to its working point, where critical oscilla- 
tions just set in, and where the hair bundle's sensitivity 
is maximal. 

The hypothesis that the self-tuning mechanism is 
based on the detection of oscillations in the transduc- 
tion current could be tested experimentally by disabling 
the electrical oscillations in the cell body. This could 
be done by voltage-clamping the cell or by specifically 
blocking, for example, the voltage-gated calcium chan- 
nels. This would disable the self-tuning mechanism and 
therefore augment the bundle oscillations. On the other 
hand, if the electrical resonance is involved in the gen- 
eration of oscillations, disabling it should eliminate the 
spontaneous bundle activity. 



Slow relaxation oscillations 

The self-tuning mechanism described above is suffi- 
ciently robust to adjust the hair bundle to the critical 
point despite the natural variation of conditions that is 
likely to occur in vivo. Many experimental 

preparations, however, involve significant perturbation 
from natural conditions, and this may shift the system 
into a different dynamical regime. In the following, we 
consider a particular type of dynamical behavior that can 
arise when the calcium feedback is weak. In this situa- 
tion, as previously discussed in Fig. EK: the interaction of 
the transduction channels with the calcium flux can pro- 
duce a bistable system. If the position xm of the adapta- 
tion motors is fixed, the hair bundle will settle at one or 
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other of the stable positions. However, this state of af- 
fairs is upset by the slow calcium feedback acting on the 
motors (Eq.EJ). Suppose, for example, that the bundle is 
at the fixed point with the higher value of x, for which 
there is a high probability that the transduction channels 
are open. The Ca 2+ ions entering through the channel 
bind to the motors, causing a fraction of them to detach, 
and the diminishing force exerted by the motors causes 
the tension in the tip links Ttl to fall. As indicated in 
Fig. Hjji, the fixed point vanishes at a critical value of 
Ttl, and the system then abruptly jumps to the other 
fixed point. At this lower value of x, the channels are 
mostly closed. The resulting drop in calcium concentra- 
tion augments the number of bound motors, increasing 
the tension in the tip links until the lower fixed point be- 
comes unstable, whereupon the hair bundle jumps back 
to its initial position. 

The cons equence of this dynamics is a slow relaxation 
oscillation l)Strogat4 II994L Sect. 7.5) of the bundle de- 
flection A, which has a non-sinusoidal form: each abrupt 
jump in position is followed by an interval of slow mo- 
tion, as shown in Fig. Eb. Martin, Mehta and Hudspeth 
(2000) have observed oscillations of this type in their 
preparations of dissected tissue. They attributed the dy- 
namical behavior to precisely the motor-mediated mech- 
anism that we discuss here, and speculated that these 
oscillations form the basis of mechanical amplification 
of auditory stimuli. According to our model however, 
the relaxation oscillations differ in many respects from 
the Hopf oscillations which provide active amplification. 
They are considerably slower; their frequency is strongly 
amplitude-dependent; and the frequency does not signif- 
icantly depend on the viscous relaxation time of the bun- 
dle, but is principally governed by the motors' sensitivity 
to changing calcium levels. In our model, relaxation os- 
cillations are an accidental consequence of the self-tuning 
mechanism when the system is too far from normal op- 
erating conditions. 

In the light of our analysis, various interpret ations of 
the experiment of Martin, Mehta and Hudspeth 120001 are 
possible. It may be that in saccular hair cells, which are 
designed to detect low frequency vibrations, slow motor- 
mediated relaxation oscillators are sufficient to provide 
a rudimentary amplification of signals. Another possi- 
bility is that the observed oscillations can be adjusted 
to the vicinity of a Hopf bifurcation by some mechanism 
that we have not considered here, e.g. by alternating 
the bundle stiffness. A third possibility to bear in mind 
is that the conditions in preparations of dissected tissue 
are not precisely the same as those in vivo, and thus 
the hair bundles may not be operating in the dynami- 
cal regime that is physiologically relevant. Careful ad- 
justment might permit the channel adaptation and the 
motor-mediated self-tuning mechanism to work together 
as we propose they do in auditory hair cells, in which case 
more rapid, critical Hopf oscillations might be observed. 




FIG. 6 Slow relaxation oscillator regime, (a) The solid line 
shows the tension-displacement relation for a channel, the 
dashed lines the contribution to tension from the passive bun- 
dle stiffness (which depends on the position of adaptation mo- 
tors). In a quasistationary state these two must be equal. For 
a given position of the adaptation motors, there can be two 
stable solutions. The motors then move the system from one 
fixed point to the other, (b) Relaxation oscillations have a 
lower frequency than and a different shape from the critical 
oscillations. Fo achieve this regime the Ca 2+ concentration 
outside the bundle was reduced and the self-tuning mech- 
anism was replaced by a constant Ca 2+ flow from the cell 
body. 



Conclusion 

Our phenomenological model illustrates that the two 
adaptation processes which have been documented in 
hair cells can, together, generate a self-tuned Hopf bi- 
furcation. The feedback mechanism that accomplishes 
the self-tuning requires, in addition, a system to detect 
bundle oscillations. We have shown that this can be 
accomplished by the resonant response of the BK and 
voltage-gated Ca 2+ channels which has been established 
experimentally in lower vertebrates. The spontaneous 
frequency u)q of a hair cell is specified principally by 
the architecture of its hair bundle. To ensure robust 
operation of a hair cell, we expect that both the cal- 
cium relaxation rate A, and the resonant frequency of 
the electrical oscillator are quite closely matched to the 
spontaneous oscillation frequency ojq. The self-tuning 
mechanism can reliably keep the hair bundle poised on 
the verge of oscillating in vivo. But the perturbations in 
conditions which can occur in some experimental prepa- 
rations may lead to quite different dynamical regimes. 
Thus, while recent experiments on hair cells of the bull- 
frog sacculus show promising signs of active amplification 
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IIMartin and Hudsoetbl ll99flL 120011: iMartin et all I2001L 

2000), they may well not yet have revealed the full po- 
tential of auditory hair cells to detect faint sounds. 



Appendix 

Series expansion 

In this section we derive an analytical expression for 
the third-order coefficient B in the complex equation of 
motion (Eq. I16|l . The Taylor expansion of Z up to third 
order, with the use of a new variable £ = Ax — Ax s (C) 
can be derived from the equations of motion (Eqs.lSJand 
EJ and the definition of Z (Ea.Hty: 

Z = [Linear terms in Ax, AC] 



iV 7 2 



C 

ccXCm 



re 



G 



(23) 



and using the specific form of Ax s 
finally obtain 



given by Eq. |3J we 



ReB = 



(A/x) 2 /iV 7 2 



16w 2 



C 



ACjvf 
C 



Mi 
D 



(1 - +D(1 + 



(A") 2 

,2 



(A M ) 2 



(27) 



If we assume a matching between the Ca 2+ relaxation 
rate and the characteristic bundle frequency, A ~ ~ 
wCj the position x of the bifurcation becomes indepen- 
dent of the characteristic frequency and so do the terms 
Po> / and all their derivatives, evaluated at the bifurca- 
tion. The third-order coefficient then scales as 



RcB- 



d 2 
"o 



(28) 



The parameters we use in our simulation give the value 
Re 6 = -&.7uj c /d 2 . 



Noise 



Expanding £ as a Taylor series in Ax and AC, expressing 
the latter in terms of Z using Eq. El 



_ ( Q + 4)Z-(q + x' s )Z 1 „ (Z - Z) 

? - - 2- 



a — a 2 s (a — a) 2 

1 ,„{Z-Zf 



6 s (a -a) 3 



+ 0{Z A ) (24) 



and collecting the real part of the terms containing Z 2 Z 
yields 



ReB = - 



7V 7 2 



16C(Ima) 2 
A Cm 



(-</"+ \a + x'ff") 



16(Ima) 



jy 4'p' -x'^2Rea + 3xM 



+ x' s \a + x' s \ 2 p'oJ . (25) 



Inserting a = x' s (— A(l+/i)+ia;o)/(A/*) from Ea.lTKl gives 



16(a. c x',) 2 ^ C 

, ,2 



(A/.) 2 



C 1 j-/// 



r +acm KVo 



- 2//x)j£ + (^) 3 (l + ^) pgA ) (26) 



In order to determine the level of the smallest stimulus 
that a hair bundle can detect, stochastic fluctuations of 
the bundle must be taken into account. There are several 
sources of noise in the system. The most important is the 
Brownian force that the surrounding fluid exerts on the 
bundle. The noise resulting from the stochastic open- 
ing and closing of channels is less significant, as can be 
seen by considering a simplified two-state channel model. 
We denote the lever arm shift between the two states by 
d, the average dwell times as To and t\ and the proba- 
bility of fining the channel in state 1 by p\. Then the 
expectation value of the number of channels in state 1 is 
(m) = Npi, and the variance is (An 2 ) = Npi(l — p\). 
The fluctuation in the total force from all channels is 
(AF 2 ) = (K T hd) 2 (Anj). The noise has the spectral 
density of a dichotomous process 



1 



N 



7T(1 + {UJTCY) 



(29) 



It has the characteristic of a Brownian noise with an ef- 

c 



fective temperature fcsTjv = ^{KT^dYNp^l— p\)rc 



ksTujcTc, where t c — (t 1 + t x 



The second ap- 



proximation is based on the estimates ksT ~ Ktl^q 
(from the condition for the negative stiffness) and luc ~ 
N~/ 2 K TL /C, (Eq. US with A ~ uj c )- Since the channel 
switching time is shorter than the period of the oscil- 
lations, ujqtc <C 1, it follows that the channel noise is 
negligible compared to the thermal noise. The same can 
be shown for the influence of stochastic channel opening 
on the Ca 2+ concentration inside the stereocilia. 

Taking the Brownian force exerted by the fluid into 
account, the equation of motion for the complex variable 
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Z becomes 



Z = (e + iujc)Z + B\Z\ 2 Z + T+^rF N (t) (30) 

where Fn (t) represents white noise with the autocorrela- 
tion function (F N (t)F N (t')) = [2(k B T/-f 2 }S{t-t'). Writ- 
ing Z' = Ze~ iuit (and taking the stimulus frequency to 
be equal to the characteristic frequency), this equation 
can be transformed to 



C -Z' = -^U{Z') + F N {t)e 



r 

U(Z' 



dZ' 

7 2 



z , 2 _ ReS / 4 _ ,j., z , 



(31) 
(32) 



It can be further simplified if we replace the "rotat- 
ing" noise FN(t)e~ luJt (which reflects the fact that the 
noise only acts on the spatial coordinate and not on 
the Ca 2+ concentration) with an averaged, homogeneous 
noise term. The averaged noise then has the tempera- 
ture |T and Eg. 1321 becomes a Langevin equation for a 
particle of friction coefficient C/7 2 > moving in a poten- 
tial U(Z'). The equilibrium solution is the Boltzmann 
distribution P(Z') oc exp(-U{Z')/\k B T). For a bundle 
that is tuned to the critical point (e = 0), the mean- 
square fluctuation amplitude in the absence of a stimulus 
{T = 0) is 



<Z 2 > 



I 2k B T 1 2 
7rC(-Re£) 



~F= ( 33 ) 



(in the second step we used the estimates Re B ~ — u> / d Q 
(Eq. I28|) and k B T ~ K^hdo, the condition for negative 
stiffness). Using the numerical parameter values we ob- 
tain 



(Z 2 ) « 0.13 



dl 



N 



Weak stimuli which evoke a response that is smaller 
than the stochastic fluctuations of the bundle, never- 
theless caus e phase-locking of the spontaneous bundle 
oscill ations ijCamalet et al. . 120001: iMartin and Hudspethl . 
l200l|) . The brain might detect such signals through peri- 
odicity in the spike-train elicited by individual hair cells, 
or via the synchronicity of spikes from neighboring hair 
cells, whose characteristic frequencies differ only slightly. 
In either case, the quantity of interest is the Fourier 
component of bundle oscillations at the sound frequency. 
With a finite stimulus force J 7 , the Fourier component is 



certain threshold value p. The smallest detectable am- 
plitude of fluid motion then scales as 



A,. 



v b 



(35) 



or with our numerical parameter values 



1.0— -TV" 
v b 



Taking p w 0.2, L/b « 3, v « 1, and N w 200 gives 
^4mm ~ 0.1 nm for high frequency bundles containing 
many stereocilia. 



Simulation of self-tuned critical oscillations 

The model parameters and their values, used in the 
simulation, are listed in Table [fl The channels were 
simulated as described in the main text, using a three- 

60 zJ, 

Onm, dc2 = 7nm, do = 8.5 nm 
and D — 4nm. The transition rate between the states 
CI and C2 and between C2 and O was modeled as 

y 



state model with parameters Gq X = 0, Gq 2 
G Q = 70 zJ, dci 



= r ch exp((Gi — Gj)/2k B T), with the free energ; 



in each state Gi — ^Ktl{x — x s (C) — d,) 2 + G° 

The resonant filter was simulated using the differential 
equation corresponding to the filter characteristics given 
by Eq. El 

V = -u 2 F (V - Vo) - ^-V + (no/N - p )Im P utRF^ 2 F 

(36) 

where no is the number of open transduction channels 
and the parameters Rf — 35mV/nA, /input = 0.5 nA, 
Q = 10, V a = -70 mV, p o = 0.15 and uo F = 2ir x 50s- 1 . 

The calcium contribution from the cell body was mod- 
eled according to Eq. (|2*2*|l with I Ca = 10^M/s, Ab = 
Is" 1 , V Ca = -65 mV, q/k B T = 0.2mV _1 . 

The adaptation motors were modeled with stochastic 
attachment at a fixed rate w m — 200 s- 1 and detach- 
ment at a Ca 2+ -dependent rate g(C)uiM, with 17(C) = 
(O.lpM- 3 ) x C 3 . 

All differential equations were solved using the Euler 
method. The time step was 1 ps. 
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TABLE I Parameter values, used in the simulation 



Tip-link stiffness 
Channel states 

Channel state energies 

Ca 2+ dependent gating point shift 
Switching time 

Distance between stereociliary roots 

Bundle height 

Geometric factor 

Pivotal stiffness of a stereocilium 

Number of stereocilia 

Viscous friction 



Ktl 
dc2 

do 

r° 
^ci 

D 

T c h 

6 
L 

7 

K 

N 

c 



0.5pN/nm 
7nm 
8.5 nm 
60 zJ 
70 zJ 
4 nm 
50 /xs 
1 fim 
7 )j,m 
0.14 

1.5 x 10 _16 Nm/rad 
50 

0.65 pN ms/nm 



(Martin ct al., 2000) 
(Martin et al.. 2000) 



(Corev and Hudspeth. 1983) 
(^jore^mclHudsDcth, 1983) 
(Howard et al.. 1988) 
(Howard et al.. 1988) 
(Howard et al.. 1988) 
(Hudspeth et al.. 2000) 
(Martin et al.. 2000) 



Ca 2+ relaxation rate 


A 


100 a' 1 


Channel contribution to the Ca 2+ concentration 


Cm 


7/jM 


Number of adaptation motors per channel 


N m 


20 


Motor association rate 




200 s" 1 


Motor speed 


v M 


0.3 (im/s 


Motor stall force 


Fm 


lpN 


Maximum transduction current 


-^Input 


/N 500 pA 


Filter resistance 


Rf 


35mV/nA 


Filter Q-factor 


Q 


10 


Filter frequency 




2tt x 50 s" 1 


Membrane potential 


V 


-70 mV 


Ca 2+ channels gating point 


V Ca 


-65 mV 


Ca 2+ channels gating charge 


9 


q/k B T = 0.2 mV 


Time constant of the slow feedback 


A B 


Is" 1 


Ca 2+ gain due to slow feedback 


I Ca2+ WfiM/s 



(Lumpkin and Hudspeth, 1998) 



(van Netten and Kros. ,2000) 
(Hudspeth and Lewis, 1988) 
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